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Abstract 

Recently an Af = 1 supersymmetric model of BCS superconductivity was proposed realizing 
spontaneous symmetry breaking of a U(1)r symmetry. Due to scalar contributions the supercon- 
ducting phase transition turned out to be first order rather than second order as in standard BCS 
theory Here we consider the effects of an external magnetic field and spatial fluctuations of the 
gap in that model. This allows us to compute the magnetic penetration length and the coher- 
ence length, and also to distinguish between type I and type II superconductors. We compare the 
supersymmetric and standard relativistic BCS results, where the main differences come from the 
different orders of the phase transition. 
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1 Introduction 

The phenomenon of superconductivity was discovered in 1911 by H. Kamerlingh Onnes, who noticed 
that electric resistivity of mercury drops abruptly to zero around 4K. Later on, W. Meissner and 
R. Oschenfeld discovered a second independent effect, now known as the Meissner effect, by which 
superconductors exhibit perfect diamagnetism. The existence of this reversible Meissner effect implies 
the existence of a critical magnetic field, H c , above which superconductivity is destroyed. 

The first description of these two phenomena was given by the London brothers with the simple 
equation J oc A for the superconducting current in terms of the electromagnetic gauge field. In 
particular, this implies that the magnetic field is exponentially suppressed inside the superconductor, 
which allows to define a magnetic penetration length, A. 

In 1950 Ginzburg and Landau (GL) proposed the following free energy 



which describes superconductivity near the critical temperature, T c , in terms of a U(l) spontaneous 
symmetry breaking. This free energy describes a second order phase transition whose order parameter 
is a complex scalar field, A, related to the local density of superconducting electrons, n s = |A(x)| 2 . 
The GL theory introduces a characteristic length, the coherence length, £, which characterizes the 
distance over which A(x) can vary without undue energy increase. Taking the quotient between these 
two characteristic lengths, n = one defines the GL parameter. 

A microscopic description of superconductivity was given by Bardeen, Cooper and Schrieffer (BCS) 
in 1957. They showed that pairs of electrons with opposite spin can bind forming a charged boson, 
known as Cooper pair, which may condense below the critical temperature. 

The same year as BCS, Abrikosov studied the GL theory in the regime k> 1, instead of k <C 1 as 
in typical classic pure superconductors. He found what he called type II superconductors to distinguish 
them from the already known k < 1 superconductors (type I). Instead of a discontinuous breakdown of 
superconductivity in a first order phase transition at H c like in type I superconductors, this new type 
of superconductor exhibits an intermediate vortex state above a lower critical field, H c ±, and below 




(1) 
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a second critical field, H C 2, above which superconductivity disappears completely. In this vortex 
state, the magnetic field penetrates in regular arrays of flux tubes of non-superconducting material 
surrounded by a superconducting current. 

In this paper we proceed with the work done in [1], in which a supersymmetric version of relativistic 
BCS superconductivity [2, 3, 4] was proposed. Apart from being an interesting problem illustrating 
the difficulties that the introduction of a chemical potential in a supersymmetric theory brings up, 
implementing the BCS mechanism within the context of a supersymmetric field theory is of great 
interest because supersymmetric theories are more stable under radiative corrections due to one- 
loop cancellations between fermionic and scalar contributions. Moreover, the fact of having scalar 
superpartners for fermions in the supersymmetric version of BCS theory can be of relevance for real 
condensed matter systems where scalar pseudoparticle excitations arise. 

The implementation of BCS superconductivity requires the presence of a Fermi surface and the 
introduction of a chemical potential for fermions, but these two conditions may lead to problems 
when scalars are included in the supersymmetric version. If the chemical potential is introduced 
through a baryonic U(1)b symmetry, supersymmetry demands that both fermions and scalars have 
the same coupling to the chemical potential. If the chemical potential is greater than the mass, a 
feature which turns out to be necessary for the existence of a Fermi surface, scalars undergo Bose- 
Einstein condensation (BEC) which prevents the BCS mechanism to be at work. Fortunately, there is 
a U{1)r symmetry which allows us to introduce the chemical potential for fermions while keeping the 
scalars neutral, avoiding in this way BEC. The main difference between this supersymmetric model 
and standard BCS, is that, as we will explain, IR physics of scalars makes the superconducting phase 
transition first order. 

Our aim here is to study the supersymmetric model in [1] under the effect of an external magnetic 
field, which will allow us to compute the magnetic penetration length and describe the Meissner effect. 
We also consider the effect of spatial fluctuations of the field A around its vacuum expectation value, 
necessary to compute the coherence length. Once known these two lengths we can compute the critical 
magnetic fields, H c \ and H C 2, that we would have for a type-II superconductor, so that compared with 
the critical magnetic field of a type I superconductor, H c , allows us to distinguish between the two 
types. 

The paper is organized as follows. In section 2, we review the relativistic BCS theory and comment 
on how to obtain the aforementioned characteristic magnitudes for a superconductor as functions of 
the temperature. In section 3, we briefly explain the supersymmetric BCS model. In section 4, we 
go on with the comparison of the magnitudes computed within the two theories, and comment on 
whether the superconductors are type I or II. Finally, we end up in section 5 with some conclusions. 
Details of the computations are given in the appendices. We will follow the notation of [7] in writing 
the Lagrangians. 

2 Relativistic BCS theory 

As explained in [4], in relativistic BCS theory one has the effective Lagrangian 

£ = ^(V'T^V - d^^ifj) - rmfjif) + [iifjty + 5 2 (V'c75^) t (^c75V') , (2) 

with a global U(l) symmetry 1 , which ensures fermion number conservation and allows the introduction 
of a chemical potential for this conserved charge. One considers the theory at finite temperature by 

1 Strictly speaking, if the symmetry is global, we are dealing with superfluidity, but transport properties are similar 
to those of superconductivity and we can consider this model as that of a superconductor in the limit in which the U(l) 
symmetry is "weakly gauged" . 
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going to Euclidean space and one performs a Hubbard-Stratonovich transformation introducing the 
auxiliary field A(x), in order to have a Lagrangian quadratic in the fields. Then the Lagrangian 
becomes 

1 i - 

Ce = -(ip^d T ip - (9 r T/>V) _ 2^*^ ~~ 9iifjj l ifj) + mipip - T 

+ 5 2 |A| 2 - 5 2 [At(v; c75 V;) + A(Vi c 75V) t ] • (3) 
From this Lagrangian we get the following fermion energy eigenvalues: 

u± = v / (co (p)±fi) 2 + ^\A\ 2 j (4) 

where ujq = \/ p 2 + m 2 , and w_ and u> + stand for particles and antiparticles respectively (supposing 
fi > 0). Once the Lagrangian is quadratic in the fermions, fermions can be integrated out to obtain 
an effective potential for A. The one- loop effective potential, V c s, which is identified with the free 
energy density in the grand canonical ensemble, is obtained by adding to the classical potential the 
thermal contribution plus a Coleman- Weinberg contribution independent of the temperature, 




written in different lines respectively. The integral over momentum for the Coleman- Weinberg contri- 
bution is divergent and one must restrict the integral to energies below the Debye one, A, which acts 
here like a UV cut-off. 

At low temperatures, the thermal contribution is dominated by the minimum of the energy eigen- 
values. This minimum is reached for w_ around ujq = /x, and contributions of uj + , i.e. those of the 
antiparticle, can be neglected. The momentum space location of this minimum defines a Fermi surface, 
pp=^ 2 — m 2 , so that n > m is required for the Fermi surface to exist. 

The condition for the formation of the fermion condensate, A 7^ 0, is that the coefficient of the 
0(A 2 ) term in the one-loop effective potential changes sign. The point where this happens is found 
by solving the gap equation, d £ V e g = 0, where we have introduced the notation e = |A| 2 . Explicitly, 
the gap equation takes the form 

o 2 f A 2 { tanh(±Pu-(p,A)) , tanh [\j3u + (p, A)) \ 

W„ pp v w_(p,a) + w+ (p,a) ; • W 

This gives the critical curve for the order parameter A as a function of the temperature, showing a 
second-order phase transition. One can obtain the value of the critical temperature by solving this 
equation at A = 0. 

If one is interested in studying electromagnetic properties, such as the Meissner effect, one must 
include in (3) a U(l) gauge field, A^. This gauge field is going to be treated as an external field. 

One can also consider fluctuations of the gap, A(x) = Ao + A(x), around its equilibrium position, 
Ao, determined by eq. (6). This will allow us to compute the coherence length, which is a measure of 
the rigidity of the condensate. For simplicity we will consider static fluctuations and suppose A to be 
real. 

We are going to treat the inclusion of the gauge field and the fluctuations of the gap as static 
perturbations. To this purpose, we write the Lagrangian (3) as C = ^O f ^, where & = O/^w), 
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and we split the matrix Of as O f (A, ^4) = Ofo(Aq) + SOf(A, A). In the path integral formalism this 
amounts to consider the saddle point approximation, which is the approach that Gor'kov [6] followed 
to derive the Ginzburg-Landau effective action from the BCS theory. In this way, the free energy can 
be expanded as 



n 



d 3 x V c i(A + A(x)) - — logdet F = + ^i + ^2 + • • • 

Zp 



(7) 



O 
Hi 



1 



d 3 xg 2 Al - ^logdetOiro , 



d 3 x2g 2 A A(x) - — Tr[Opl50 F ] 
n 2 = J d 3 xg 2 A 2 (x) + ^Tt[(O f 150f) 2 } 



d 4 x J d 4 xitr[OQ 1 (x,x 1 )50(xi,x)} , (8) 
d 4 x / d 4 xi / d 4 x 2 / d A x 3 tr[OQ 1 (x,x 1 )50(xi,X2)OQ 1 (x2,x 3 )50(x 3 ,x)] , (9) 



where 

Tr[O^SO] = 
TrKO^dO) 2 ] = 

50(x,y) = 5^(x — y)dO(x) and tr[.] is the usual trace over matrix elements. The first term of the 
expansion (7), Oo, is just the free energy corresponding to eq. (5), which fixes the value of the gap. The 
second term, Qi, has two contributions, one corresponding to A (which vanishes, as it is proportional 
to the gap equation) and another one due to the gauge field A. This last contribution gives rise, in 
momentum space, to terms of the form A ■ p, which vanish by isotropy once integrated. This just 
leaves contributions involving the temporal component of the gauge field, which are interpreted as 
fluctuations or space inhomogeneities of the chemical potential. 

Once the effective potential is computed we have to add the kinetic term for the gauge field, the 
complete free energy is then 



n 



tot 



d 3 x 
d 3 x 



-E 2 + -B 2 
2 2 



+ n(A,A) 



fl 



V - ieA I A 



+ m- 2 \A\ 2 + f 2 E 2 + f 3 B 2 



+ tt + ... 



(10) 



where the coefficients f±, f 2 , f 3 and m~ 2 have been introduced to account for the contributions coming 
from the Oi and ^2 terms of the expansion of the free energy. Near T c , we can expand 0,q in eq. (10) 
as 



tto(A ,T) 



d 3 x 



a(T - T C )A 2 + bA A + ... 



(11) 



where 



^~2^c J d P P 2 ^sech 2 Q/3 c (w (p) + fJ>)) + (ji -> , 

sech 2 (±/3 c (u (p) + n)) tanh (\p c (u) Q (jp) + //)) 



(u (p) + /x) 2 



+ 2- 



(uj (j)) +/x) 3 



+ (a* -»• -m) 



This is just the GL free energy, which is valid near the critical temperature, so (10) must be considered 
as a generalization of the GL free energy, valid for a wider range of temperatures. 



5 



From the equations of motion for A (in the London limit, i.e. A = 0, and in the Coulomb gauge, 
V • A = 0) one sees that the magnetic field is exponentially suppressed inside the superconductor, 
which is the footprint of the Meissner effect. The corresponding magnetic penetration length A is 
given by 

V 2 A = ^\A \ 2 A => ^ = £^l|Ao| 2 . (12) 

J3 * J3 

Another characteristic length is the coherence length, £, which characterizes the distance over 
which the superconducting electron concentration cannot drastically change. This length is obtained 
from the equations of motion for A(x) obtained from (10). In absence of gauge field, we get 

Hence to determine these magnitudes we must be able to identify the coefficients fi, and m~ 2 
in (10). The identification of these coefficients is explained in appendix B. 

Other magnitudes we want to compute are the specific heat and the critical magnetic field above 
which superconductivity is destroyed. The former was already computed in [1] at constant /i, instead 
of constant charge density, as it is a more sensible quantity when scalars are included in the SUSY 
case. The specific heat is computed through the formula 

_ fdV cS \ r r fdS\ fd 2 V cS d 2 V cS de\ de drdeVeS nA , 

S = -Vdr)/ c = T \lr)=- T \-W^ + Wd-e&T)' where df = ' (14) 

The critical magnetic field is obtained by equating the work done for holding the magnetic field out 
of the superconductor with the condensation energy, 

&l = V n (T)-V s (T) , (15) 

where V n and V s are the free energies per unit volume in the normal and superconducting phase at 
zero field. This is the critical magnetic field over which superconductivity would be destroyed in a 
type I superconductor. For a type II superconductor, there are two critical magnetic fields, in between 
we find the vortex state. H c \ is the value of the magnetic field for which a single vortex with flux 
quantum <j)Q = 1T /e would appear. Whereas near H C 2 vortices are as closely packaged as the coherence 
length allows. Therefore, 

Ha « . H *~^e- (16) 



3 Supersymmetric BCS theory 

The supersymmetric version of BCS theory [1], is based on an M = 1 SUSY theory described by the 
following Kahler potential, 

K = &$ + g 2 (&$) 2 , (17) 

and no superpotential. In [1] it was necessary up to four chiral superfields. Two superfields with 
the mentioned Kahler potential for each one, in order to have Dirac fermions, just as in standard 
BCS theory; and two more free superfields with canonical Kahler potential to ensure that the U{\)r 
symmetry, for which the chemical potential is introduced, is non- anomalous. As these four chiral 
superfields are decoupled form each other, it is sufficient to study the one superfield model, (17). 

The chemical potential is introduced in this theory for a U(1)r symmetry, in such a way that 
scalars are neutral and, in consequence, fermions have .R-charge —1. This is done in this way in order 
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to prevent the scalar from coupling to the chemical potential, thus avoiding Bose-Einstein condensation 
(BEC) of the scalar. In this minimal setup, with one chiral superfield, this particular choice of R- 
charges prevents the introduction of a superpotential, which may lead to mass terms, in this way we 
are dealing with a massless theory. 

This theory also has a baryonic U(1)b symmetry. If the chemical potential is introduced for 
this U(1)b symmetry instead for the U(1)r symmetry, this would inevitably lead to BEC for the 
scalar in the presence of a Fermi surface. In this case, there is no reason why we cannot introduce 
a superpotential that generates mass terms, for example W = m&&. With this superpotential we 
obtain the following scalar energy eigenvalues, 



lo s ± = Vp 2 + m 2 + 4# 4 A 2 ± H ■ (18) 

As we have seen, the existence of a Fermi surface requires [i > m, for which ujs- would become negative 
at low momentum near the non-superconducting phase (A = 0), making the thermal contribution to 
the effective potential ill-defined, resulting in the appearance of BEC. 

Once expanded the Kahler potential (17) in terms of component fields, we have to perform, as 
in the previous section, the Hubbard-Stratonovich transformation introducing the auxiliary field A, 
which now is related with the fermion condensate through 

* = r-rW- ( 19 ) 

1 + Ag 2 v 2 

After expanding the scalar field to quadratic order around its vacuum expectation value, v, which turns 
out to be zero, going to Euclidean space and introducing the chemical potential for the i?-symmetry 
we arrive at the following Lagrangian 

C E = (d a (f)d a 0* + {d T + /i) V> + i^a'Sty) + 4 5 4 |A| 2 |0| 2 - g 2 A(W) - g 2 A*$$) , (20) 

with classical potential V c \ = g 2 |A| 2 . 

Now we want to study the magnetic response of the superconductor, therefore we have to turn 
on an external U{\) gauge field. We have two possibilities: the gauge field can be turned on for the 
baryonic U(1)b symmetry or the U(1)r symmetry 2 . In either case supersymmetry is broken by the 
background. In principle, the magnetic response could depend on this choice, but as explained at the 
end of appendix B, this is not the case, at least for a small enough gauge coupling. For definiteness 
we introduce the gauge field through the baryonic U(1)b symmetry, therefore the Lagrangian (20) is 
modified to 

C E = (d a (f)d a (j)* + i)a° {d T - ieA T + p) y> + i$d* (d i - ieA 1 ) if>) 

- ieA a 4>d a (j)* + ieA a <j)*d a (t> + e 2 A a A a \(j)\ 2 + 4 5 4 |A| 2 |0| 2 - g 2 A(^) - g 2 A*(^) . (21) 

Given this Lagrangian, we can construct the bosonic and fermionic matrices (whose explicit form is 
written in appendix A) and split them as explained in the previous section, so that we have to add to 
eq. (7), the scalar contribution 

+ -L logdet O S0 + ^Tr[O^SO s ] - ^Tr[(0^80 s ) 2 ] + ... (22) 



2 -R-symmetry can only be gauged within the context of supergravity. The present model can be easily embedded in 
M = 1 supergravity, in such a way that in the usual laboratory set up, where energy configurations are much lower than 
the Planck scale, the supergravity multiplet can be ignored. 
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The determination of the magnitudes explained in the relativistic BCS theory is now completely 
analogous for the supersymmetric case. For example, the effective potential analogous to (5), becomes 
now 

V cS = g 2 A 2 



1 r d 3 p 

~PJ (2^)3 
1 f d 3 p 



2 J (2tt) 
where the energy eigenvalues are 



/ (2^3 ( bg(1 + e ~^ +) + l0g(1 + e "^ _) " 21 ° g(1 " e ~^ S) ) 

3 (w + + w_ - 2lu s ) , (23) 



oj f ± = ^(p±^Y + 4g*A 2 , oj s i,2 = \jv 2 + VA§ . (24) 

4 Comparison between the SUSY model and relativistic BCS 
4.1 Gap 

According to the previous section, we have to solve the gap equation d £ V e R = 0. Just as shown in [1], 
the explicit form of the gap equation is 



g 2 ,A /tanh(i/V 4 9 4 Ag + (p-M) 2 ) tanh(j/V4<7 4 Ag + (p + /i) 2 



2coth(±/VVA ; +p 2 



VVAgTp 2 



(25) 



where the second line is the new contribution due to the scalar. For the standard relativistic BCS 
case the effective potential develops a non trivial minimum at Ao 7^ 0, which will approach the origin 
as we increase the temperature reaching the zero gap value at the critical temperature T c , (fig. 1 (a) 
and (b)). The phase transition is second order. For the SUSY case, at low temperatures, T < T c \, 
the effective potential has a unique non trivial minimum (zone 1 of fig. 1 (c) and (d)), above T& (zone 
2) a zero gap metastable minimum appears becoming the dominant one for T > T C 2 (zone 3), so that 
the minimum with non vanishing gap becomes metastable until it disappears at T C 3, jumping to the 
non-superconducting zero gap solution for T > T C 3 (zone 4). Therefore, the phase transition is first 
order in the supersymmetric case. For short, we will be calling T c = T C 3 and A c = A(T C 3) from now 
on. 

As explained in [1] , the IR physics of massless scalars make the phase transition first order rather 
than second order. First order phase transitions in usual superconductivity have already been ex- 
plained in [5], where a gauge field takes a non zero vacuum expectation value appearing as a ( J 4 2 )|Ao| 3 
term in the free energy, which inevitably leads to a first order phase transition. This gauge field plays 
the same role as the scalar in the supersymmetric case. Near Ao = and at low momentum we can 
approximate the contribution of the scalar as 

2g 4 T /■_,_ 2 1 2g e T 



deV eS ps — — / dpp Pi A , (26) 

which leads to the analogous C(|Ao| 3 ) term in the free energy. 

At low temperatures, we can neglect in (25) the antiparticle and the scalar contributions for 
momenta near the Fermi surface, where the main contribution to the integral comes from. By doing 
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Figure 1: (a) Effective potential as function of the gap for different temperatures in the relativistic 
BCS case, (b) Corresponding gap as function of the temperature, (/i = 0.05, g = 0.07, A = 65). (c) 
and (d) analogous figures for the SUSY case, (p, = 5, g = 0.5, A = 400). The thicker red lines in (d) 
correspond to metastable solutions and the curved dashed line represents the maximum that separates 
both minimums in the effective potential. 



so, we are taking the non-relativistic limit and connecting with the standard BCS result. Specifically, 
one performs the ad hoc approximation by which one substitutes dp 3 by A-Kp 2 F dp and the integral is 
done in the interval \p — pp\ < A around the Fermi surface. After these approximations, one can show 
that the behaviour of the gap near the critical temperature T c follows a universal behaviour 



where r\ = 1.74 can be computed numerically from eq. (11). In the supersymmetric case we obtain 
numerically the following behaviour near the critical temperature 

(28, 



Ao(0) ' V T Cj 

where a ~ 0.5 and the universal behaviour is lost, since now rj depends on the parameters g, \x and A. 

One of the virtues of supersymmetry is the softening of divergences due to cancellations between 
fermionic and bosonic contributions. This feature is also observed for the present system. As the ther- 
mal integrals in (5) and (23) are convergent, it is sufficient to study the dependence of the gap with the 
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Figure 2: Specific heat a function of the temperature, (a) Relativistic BCS. (b) SUSY BCS. 

cut-off at zero temperature. The dependence in both the relativistic BCS theory and supersymmetric 
BCS theory is 

A 3 A -- s^-S 



|A (0) re l.BCs| ~ ^2 



o(0)siJSYBCS| 



(29) 



where we have assumed A S> <? 2 Ao, /i. Thus we see that the cut-off dependence is milder in the SUSY 
case. If one performs the aforementioned approximation to consider the non-relativistic limit, the 
dependence changes to 



|Ao(0) 



A 



BCS I 



2tt 2 
7 — 1 

e 9 /* 



(30) 



The dependence in the non-relativistic case is linear rather than cubic as in the relativistic case. This 
reduction in the power of the cut-off is due to the substitution of dp 3 by A-Kp 2 F dp. 



4.2 Specific heat 

The comparison between the SUSY and the relativistic BCS case for the specific heat, computed with 
eq. (14), is shown in fig. 2. When only fermions are considered the jump in the specific heat at the 
critical temperature is finite, as it is characteristic for second order phase transitions, whereas it is 
infinite when the scalar is included. 

The behaviour of the specific heat in the different temperature regimes is explained in [1], whose 
expressions we show here 3 . At high temperatures (Ao = 0) the different contributions to the specific 
heat corresponding to the scalar, the fermion and its anti-particle, are the following ones 



cs|a =o 



4vr 2 T 3 
15 



cf±|a =o 



7vr 2 T 3 
60 



for T » fj, . 



(31) 



At low temperatures the antiparticle and the scalar can be neglected and the main contribution 
comes from the region near the Fermi surface p ~ /x, so that the behaviour of the specific heat in the 
normal phase and in the superconducting phase is given by 



cf-|a =o 



ji 2 T 
6 



cf-|Ao^o ~ e 



2<^A (T=0) 
T 



(32) 



3 There is a 1 /2 factor between the results presented here from those of [1], because here we are considering just one 
complex scalar field and a Weyl fermion. 
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The last expression shows a way to compute the value of the gap at zero temperature. Indeed, by 
performing a fit of the plots in fig. 2 one can obtain the value of the gap at zero temperature shown 
in fig. 1. 



4.3 Magnetic penetration length and coherence length 

The magnetic penetration length, obtained from eq. (12), is plotted in fig. 3. In both the relativistic 
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Figure 3: Magnetic penetration length as a function of the temperature, (a) Relativistic BCS. 
(b) SUSY BCS. 

and the SUSY BCS theory, the magnetic penetration length is a monotonically increasing function 4 . 

For the fermionic contribution alone, the magnetic penetration length diverges near the critical 
temperature. Whereas in the supersymmetric case, it reaches a finite value at T c and it jumps to 
infinity for T > T c , which is the expected behaviour for a first order phase transition. This behaviour 
is explained basically from the gap dependence. Expanding the /i coefficient as a power series of the 
temperature and the gap, one obtains 



fi rd.BCS = a + oiAq + a 2 (T -T c ) + ... 



fl SUSY bcs = «o + «i(Ao - A^) + 



(33) 



where one can check that the coefficients shown do not vanish. The ellipsis stands for higher powers of 
the temperature, taking into account the temperature dependence of the gap, (27) or (28). Substituting 
this expansion in (12), and using (27) or (28), one finds the following behaviour 



Kci. 



BCS 



T 



-1/2 



(A - Ac) 



SUSY BCS 



T 



(34) 



The behaviour of the magnetic penetration length at zero temperature can be computed analyt- 
ically. According to the dependence of the coefficients f\ and m~ 2 with the Debye energy, given in 
appendix B, and that of the gap, we find the following expressions for the magnetic penetration length, 



A 



rcl. BCS 



4g 3 



3x/37r 2 



-A 



Asusy 



2tt 



bcs 



ec 



(l+4c?) 3 / 4 A -1 , where c 



exp 



7T 



g 2 ^ 2 



. (35) 



The coherence length, £, obtained from eq. (13), is plotted in fig. 4. It is a monotonically increasing 
function of the temperature for both the relativistic and the SUSY BCS theory. The behaviour near 



4 For high enough values of the chemical potential and the cut-off a counterintuitive non-trivial minimum can appear 
before reaching the critical temperature, which would mean that there is a range of temperatures where the Meissner 
effect is enhanced with increasing temperature. This odd behaviour can be avoided by restricting the parameter range 
of validity to not very high values of the chemical potential and the cut-off. However, this restriction is relaxed in the 
SUSY case where the cut-off dependence is softened. 
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the critical temperature is the same as for the magnetic penetration length. To see this, expand the 



m 



coefficient as we did with the f\ coefficient. Using the gap equation, one can see that the m 



coefficient has a global Aq factor so that the expansions are 



in 



rel 



; BCS = A 2 (b + b 1 A 2 + b 2 (T-T c ) + . 



m SUSYBCS = A§(A> + /3i(A§- A*) + . 



(36) 



Inserting the expansions (33) and (36) in the expression for the coherence length, (13), we find the 
following behaviour 



£rel. 



BCS 



1 



T 



-1/2 



(£-6 



c^SUSYBCS 



1 



T 



(37) 



Thus, the coherence length exhibits the same behaviour as the magnetic penetration length. However, 
the behaviour of the coherence length at zero temperature is different from that of the magnetic 
penetration length. Now, the dependence with the Debye energy is 



£susy 



BCS 



6 



rel. BCS 



9V6V 



A" 



2(1 + 4c 2 ) 3 / 4 



o 2 , 22 32c 4 + 16c 2 + 5 l + ^/T+4? 

27r + 9 V T-, , : o^/o 2 log - 



-1/2 



(1 + 4c 2 ) 5 / 2 



2c 



+ 0(A~ 



(38) 



(39) 
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Figure 4: Coherence length as a function of the temperature, (a) Relativistic BCS. (b) SUSY BCS. 



If we take the quotient between these two characteristic lengths, 



A _ 1 
K ~ 1 ~ v^em/iAo 



(40) 



we get the Ginzburg-Landau parameter, shown in fig. 5. As the coherence length behaves in the same 
way as the magnetic penetration length near the phase transition, at leading order, the GL parameter 
will take a finite constant value. Depending on the value of the GL parameter one has a type I (k C 1) 
or a type II superconductor (k > 1). In the GL theory k is defined near the critical temperature and 
the critical value differentiating between the two types of superconductor is k = 1/ y/2. As shown in 
fig. 5, k > 1, since we are considering small values of gauge coupling e, then the superconductors are 
type II in both the relativistic BCS case and the SUSY case. 
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Figure 5: Ginzburg-Landau parameter as a function of temperature, (a) Relativistic BCS. (b) SUSY 
BCS. 



4.4 Critical magnetic fields 

The critical magnetic field is obtained by equating the energy per unit volume, associated with holding 
the field out against the magnetic pressure, with the condensation energy. That is eq. (15), 



H 2 C (T) 
8vr 



V n (T)-V s (T) . 



The behaviour near the critical temperature is found by preforming expansions similar to those 
for the fi and m~ 2 coefficients. In the standard relativistic BCS theory the gap is expanded around 
Ao = 0. On the other hand, when the scalar is considered, H c does not make sense above T C 2, since 
the superconducting minimum in the effective potential becomes metastable, but we can perform the 
expansion around A C 2 = Ao(T C 2). According to the dependence of the gap with the temperature (27) 
in the relativistic BCS theory and due to the fact that Ao is linear with the temperature near T c2 in 
the SUSY BCS theory, we have 



u 2 

n c reL BCS 

8vr 



d T d £ {V n (T)-V s {T)) 



T=T (T-T C )A 2 



+ l -d 2 £ {V n {T)-V s {T)) 



A 4 + 

T=Tc t • • • 
A =0 



SUSY BCS 
8vr 



dr{V n {T) - V S (T)) 



T=T c2 
A =A c2 



(T - T c2 ) + . . . 



(41) 
(42) 



from which we find a linear and square root behaviour near the critical temperatures, 



H, 



c rel. BCS 



1 



T 



H, 



c SUSY BCS 



1 



T 



(43) 



as shown in fig. 6 (a) and (b), respectively. The expressions for the critical magnetic field at zero 
temperature are 



H, 



c rel. BCS 



2t5 A 3 



3tt 2 



H, 



c SUSY BCS 



g 2 fi 2 



1 



. 1 + y/1 + 4c 2 \ c . 

_ — log -A . 

n 2 \^l + Ac 2 2c jg 

(44) 



As expected the dependence with the Debye energy is milder in the super symmetric case. 
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As explained at the end of sec. 2, for type II superconductors, which are characterized by the ap- 
pearance of Abrikosov vortices in a mixed superconductor-normal state, there are two critical magnetic 
fields, H c \ and H C 2, 

<t>o 



H, 



cl 



H, 



2vrA 2 ' " c2 ~2< 2 - 

They are plotted in fig. 6 together with H c . It is easy to obtain the behaviour in the different regimes 
of these two critical magnetic fields using the expressions for the magnetic penetration length and the 
coherence length. 

If H c <C H C 2, we will be able to see the intermediate vortex state as we decrease the applied 
magnetic field, i.e. the superconductor is type II. On the contrary if H c S> H C 2, we reach the pure 
superconducting state without the formation of any vortex, and we have type I superconductivity. 
According to fig. 6, we have type II superconductivity in both the relativistic and the SUSY BCS 
theory, in agreement with the prediction obtained in the previous section by computing the GL 
parameter. 

Given that H c ends at T C 2 instead of T C 3, as H c \ and H C 2, and that there is a range of parameters 
where H C 2 < H c at zero temperature, one can find a crossing between the two magnetic fields and a 
crossover between type I and type II behaviour as we increase the temperature. However, this crossing 
effect disappears if the gauge coupling is sufficiently small. 
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Figure 6: Critical magnetic fields H c (black), H c \ (blue) and H C 2 (red) as a function of the temperature, 
(a) Relativistic BCS (e = 0.2), for clarity H c \ and H C 2 have been rescaled by a factor 500 and 0.02 
respectively), (b) SUSY BCS (e = 0.2). 



5 Conclusions 

In this paper we have computed some relevant magnitudes in superconductivity as a function of the 
temperature for the super symmetric BCS model proposed in [1] and for standard relativistic BCS 
theory. These magnitudes are the gap, the specific heat, the magnetic penetration length, which 
describes the Meissner effect, the coherence length and the different critical magnetic fields. The 
results obtained in both theories have been compared to clarify the role of super symmetry. To compute 
these magnitudes we have considered spatial fluctuations of the gap and introduced an external gauge 
field for the baryonic U(1)b symmetry. Both the gap fluctuations and the gauge field are introduced 
as perturbations, so once the effective potential is expanded and the coefficients in the generalized 
Ginzburg-Landau equation (10) are identified, the computation of the aforementioned magnitudes is 
straightforward. 
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Considering the quotient between the magnetic penetration length and the coherence length, which 
defines the GL parameter, shows that the studied superconductors are type II. This is confirmed when 
studying the relation between the critical magnetic fields H c and H C 2. 

The main differences between the standard relativistic BCS results and those from the supersym- 
metric model lies in the different orders of the phase transition. For example, the infinite jumps found 
in the supersymmetric model for the specific heat, magnetic penetration length and coherence length 
are those characteristic of a first order phase transition. Another important difference is the softening 
in the cut-off dependence in all quantities thanks to super symmetry. 
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A Scalar and Fermionic matrices 

Writing the quadratic part of the Lagrangian in the presence of a baryonic gauge field, (21), as 



C = &O s $ + &0 F V , with $ 



Pi-P) 



and ^ 



/ Mp) \ 

4(-p) 
V 4(-p) J 



the scalar and fermionic matrices after the splitting 0(A,A) = Oq(Aq) + SO(A,A), are given in 
momentum space by 



Oso 



1 (, .2 



+ p 2 ) + 2g 4 A 2 







\{u 2 +p 2 )+2g i A 2 I ' 



SOs 



( e(A T uj+A-p) + ±e 2 (a 2 . + A 2 
+2g 4 (2AA + A 2 ) 





FQ 



50 F 



V 




-5 2 A 







\u-\p+\n g 2 A 



-e(A T io + A • p) + \e 2 [A 2 T + A 2 
+2g 4 (2AA + A 2 ) 
-g 2 A \ 



V 






-g 2 A 



9 2 A 




5 2 A 



l u} + h- \n 



2 ^ 1 2 y 2 

jw-^-^y 

^2 







-5 2 A \ 



1 



(45) 



\eA T - ±eA J 

If instead of introducing a baryonic U(1)b gauge field, we had introduced U(1)r gauge field, the 
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scalar and fermionic matrices would have been the following ones: 



JOs= (V(2AA + 



50 F 



A 2 ) 

2g 4 (2AA + A 2 

/ \eA T + \eA -g 2 A 

' \eA T -_\eA g 2 A 





V 



-ff 2 A 



5 A 




7 eA T 











(46) 



The energy eigenvalues computed for the O50 an d O^o matrices are 



u F ± = J( P ± t i) 2 + 4g*A 2 



W 5 1,2 = \/p 2 + % 4 A 2 . 



(47) 



B m 2 , /1 and ^3 coefficients 

In momentum space, the m -2 and /1 terms in (10) are given by 

d 3 k 



J 



d 6 xm~ z AA* 



(2^)3 
(2tt) 3 



m- 2 A*(fc)A(fe) , 



hk'A*(k)A(k) , 



(48) 
(49) 



where we have considered time independent perturbations. Thus, we have to find in (7) a term 
quadratic in A and expand its coefficient up to quadratic order in momentum. The zero order term 
will correspond to to" 2 and the coefficient of the quadratic term in momentum will be identified with 
fi. Terms quadratic in A are found in Q 2 and, if scalars are considered, in Q^. The term only 
contributes to the to -2 coefficient, which in momentum space becomes 



1 



1 



2 ^[O^Os]- 2p2j ~^j (27r)3 



d 3 K 



tr 



0^(u n ,K)SO s (x) 



(50) 



where we have to sum over Matsubara frequencies, u n = 2mr/{3 for bosonic frequencies and co n = 
(2n + l)7r//3 for fermionic ones. Once the Matsubara sums are done, we have to consider the piece 
quadratic in A (supposing A to be real) 



If f d 3 K 
-Tr W^Os\ KK = ^5a®AA(*)A(S)A 



(x) 



1 f d 3 K 

Hence, we identify the first contribution to to -2 as 

n 1 f d 3 K 



d 3 k 
(2^)2 



A*(k)A(k) 



TO 



PJ (2tt)5 



(51) 



(52) 



Let us elaborate now on the ^2 contribution, 



^Tr[(0^50 F ) 2 ] = ^J d 4 Xl J d 4 x 2 tr[dO F (x 1 )O F 1 (x 1 ,x 2 )5O F (x 2 )O F 1 (x 2 ,x 1 )] , (53) 
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plus the analogous scalar term if one considers the supersymmetric case, and the extra term f d 3 x g 2 A 2 
for m~ 2 . The previous expression in momentum space is 

d 3 ki d 3 k 2 d 3 qi d 3 q 2 



(2tt) 3 (2tt) 3 (2tt) 3 (2vr) 3 

x e -i(u m -u n )(T 1 -T 2 ) e -i(ki-k2+$i)3i e -i(-k 1 +k2+q2)x2 

xti[SO F (qi) Op] {0J m M)5O F {q 2 ) Opl {uj n ,k 2 )\ 
tv[50 F (k 2 - h)Opl(u n ,ki)80 F (ki - k 2 )0 F l(uj n ,k 2 )} 



l ^ f d 3 k x d 3 k 2 
4/3 2-> J (2^) 3 (2^) ILi 



l f d 3 h d 3 k 2 -> -> 

F[ki,k 2 ) . (54) 



^Iaa = S I 7^7^ A *( fc 2 " h)A(k 2 - fci).F AA (fci, fe) . (55) 



^| AA = ~ a I , ^ 3 ^ 3 A^)A(fe)7 M (^,fc) , (56) 



?7t 2 



~ /3 7 (2vr) 3 (2tt) 3 

Taking the piece quadratic in A in (54), we have 

l f d 3 h d 3 k 2 
pj (2vr) 3 (2tt) 3 

Assuming that second order corrections are located close to each other in momentum space, we can 
expand the momenta around their average value, k\ = K — k /2, k 2 = K + fc /2, so that the corresponding 
free energy term is 

l r d 3 K d 3 k 
pj (2vr) 3 (2vr) 3 

Expanding U P to quadratic order in k, we identify the m~ 2 and f\ coefficients with 

l r d 3 K - 
9 2 + p J t^3(^aa(K,0) - B AA (K,0) + <B^(K)) , (57) 

1 f d 3 K 

h = WJ (^3(^a(^0) - d 2 B AA (K,0)) , (58) 
once the analogous bosonic contribution is included. 

The /3-term in the Ginzburg-Landau free energy (10) will have a contribution coming from the 
gauge field kinetic term plus contributions coming from the part of £l 2 quadratic in the gauge field, 
which will be proportional to the square of the gauge coupling, e 2 , 

h = l - + 0(e 2 ). (59) 

As the gauge coupling is assumed to be small we can simply take f% = 1 /2. From this identification 
for the /3 coefficient we see that there is no significant difference between turning on a baryonic gauge 
field (45) or an i?-symmetry (46) gauge field, since differences would appear to order 0(e 2 ). 

Once these coefficients are computed, one can study their cut-off dependence in both the relativistic 
and supersymmetric case. We will restrict ourselves to the zero temperature regime, where integrals 
can be performed analytically. In the zero temperature limit, the explicit forms of the coefficients f± 
and m~ 2 are 

128<? 12 A<j^ + 80g 8 A^(A + V) 2 + A^( A + /") 5 + 2 5 4A o( A + A 4 ) 3 (3A 2 + 6A/x + 8/j 2 ) 



fi 



T= ° 967r2Ao(4 5 4A 2 + (A + ^)2) 5/2 

+ (A* -> ~a0 

<7 8 A 2 A 3 



2tt 2 (4 5 4 A 2 + A 2 ) 



5/2 ' 



(60) 
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m 



t 
4vr 2 



I 2| g 4 / (A + /i)(5^ + 2A/i-A 2 ) + 4g 4 (5^-3A)A 2 

lr=o 51 + 



+ 2 Ou 2 - 6/ A 2 ) log ^ 

g 4 / 12g 4 A 2 A + A 3 
2^ ^y^A^ + A 2 



V4 5 4 A 2 + (A + ^)2 
^ + v / 4g 4 A 2 + j » 2 



+ 



A + /i + V4g 4 Ag + (A + A t)2 < 
-12^A 2 csch-(^)) 



+ -> -/i) 



(61) 



where the last line in (60) or (61) corresponds to the scalar contribution. Taking into account the 
cut-off dependence of the gap at zero temperature, (29), we find the following expressions for the j\ 
and m -2 coefficients at leading order in A: 



1 rel. BCS 



3V 



A 



-10 



/ 



1 SUSY BCS 



;(l + 4c 2 r 3 / 2 + d?(A- 



™J.bcs = 9 2 + 0(A- 



m 



SUSY BCS — » 1 2tt 2 
2, 



g 4 fi 2 ( 3 2c 4 + 16c 2 + 5 

(1 + 4c 2 ) 5 /2 



2 log 



1 + Vl + 4c 2 
2c 



+ 0{A~ 



(62) 



(63) 



where c = exp —-^zjp — f • We must stress that the coefficient f\ in the relativistic BCS theory does 
not vanish at zero temperature, because A is a physical cut-off, the Debye energy, which takes a finite 
value. 
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